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Abstract 

We give numerical evidence that the location of the first order phase 
transition between the low and the high density phases of the one dimen- 
sional asymmetric simple exclusion process with open boundaries becomes 
sample dependent when quenched disorder is introduced for the hopping 
rates. 
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The effect of quenched disorder on phase transitions is one of the most 
studied and best understood aspects of the theory of disordered systems 
PEIEIIIISIEIITIIHIEI- For equihbrium systems, the Harris criterion gUS] 
allows one to decide whether the critical behavior of a second order phase 
transition is altered by a weak disorder. One also knows [3 El El that 
first order transitions are suppressed by a weak disorder in low enough 
dimension. 

It is now well established that non equilibrium systems can exhibit phase 
transitions even in one dimension (for a review see ^01)- For simple systems 
as the asymmetric exclusion with open boundaries (the TASEP) , the phase 
diagram is known exactly [111 I12L 1131 114L [T5| . It is then a natural question 
to investigate what is the effect of disorder on such a phase diagram [TBI 15]. 

The one dimensional asymmetric simple exclusion process fASEPl|17[ 
I18lll9j has been considered in several contexts ranging from biopolymerisa- 
tion and electrophoresis to the study of car traffic [201 ■ I* is also related to 
other problems of statistical physics like surface growth and directed poly- 
mers in random media |21[ . It describes a one dimensional lattice of L 
sites, each site of which is either empty or occupied by a single particle (ex- 
clusion rule). During each infinitesimal time interval dt, a particle at site i 
attempts to move to its right neighboring site with probability pdt and the 
move is successful only if the target site i + 1 is empty. At the boundary 
(site 1 and site L) the dynamics is modified as follows (open boundary 
conditions): during each infinitesimal time interval dt, a particle is injected 
on site 1 of the lattice with probability adt if this site is empty; and if a 
particle is present on site L, it is removed with probability (3dt. When the 
boundary parameters a and /3 are varied, the model exhibits several phase 
transitions m^l^SElEl^ second order phase transitions along the 
lines a = p/2 < (3 (from a low density phase to a maximal current phase) 
and [3 = p/2 < a (from a high density phase to a maximal current phase), 
and a first order transition from the low to the high density phase along 
the line 

a = /3<| (1) 

The goal of the present work is to study the effect of quenched disorder 
on this first order phase transition. Different ways of introducing quenched 
disorder have been considered for the ASEP. The most studied is the parti- 
clewise disorder |23lll6ll^ . It corresponds to the case where the hopping 
rates p depend on the particle which attempts to jump. One can think 
of it as a model of traffic along a one lane road where each vehicle has 
a different, random velocity. This model can be exactly solved through a 
mapping to the zero-range process . 

We consider here the less studied sitewise disorder, where to each lat- 
tice site is associated a random, quenched hopping rate pi. In the traffic 
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analogy, this would correspond to the case where localized accidents or bot- 
tlenecks locally decrease the mean speed of vehicles. Even the adjunction 
of a single localized defect in the hopping rate Pi transforms the ASEP into 
a difficult, and yet unsolved problem |25l I26L I77| . There are only few an- 
alytical results in the case of site disorder, the most noticeable ones being 
the existence of an hydrodynamic limit for the density, the concavity of the 
current of particles as a function of the density of particles j2H| and the 
reflection invariance of the current [221 ■ Monte-Carlo simulations on ASEP 
on a ring geometry (i.e. with periodic boundary conditions') |3UI l?T| have 
shown that the current-density relation is modified by quenched sitewise 
disorder, with the appearance of a plateau between two densities Pc and 
1 — Pc where the current becomes independent of the density of particle p 
on the ring. Tripathy and Barma |30| have also shown that this plateau 
corresponds to a segregated- density regime where the disorder induces a 
phase separation between a high density and a low density phase. Outside 
this plateau (i.e. for densities p < pc or p > 1 — pc), the steady state profile 
is uniform on a macroscopic scale (but with microscopic variations). Some 
bounds for pc have been obtained jJHj. When the distribution of the pi 
allows values as close as possible to a third regime occurs, in which the 
steady state current vanishes in the thermodynamic limit ^1 OOl • 

In the open boundaries case, when particles are injected at site 1 and 
removed at site L on the last one, one expects the system to have the same 
3 phases (low density, high density and maximal current) as in the pure 
case. The question then is how the nature and the location of the phase 
transitions are modified by the effect of disorder. 

Here we investigate numerically the effect of disorder on the first order 
transition (QJ. We choose a binary distribution of the bulk hopping rates 



with Pniin = 0.8, Pniax = 1.2. 

In equilibrium systems (with short range inteactions), the properties 
of finite systems depend smoothly on external parameters and phase tran- 
sitions usually occur in the thermodynamic limit, i.e. in the limit of an 
infinite system size. The ASEP on a finite chain with or without disorder 
is a Markov process with a stationnary state which also depends smoothly 
on the parameters a and P and one needs to take the limit of an infinite 
size to observe phase transitions. 

As the Pi's depend on i, one needs to choose a procedure to add new 
sites to a given sample in order to take the thermodynamic limit. In the 
following we compare two different procedures: 



Pi- 




Pmin with probability ^ 
Pmax with probability ^ 



(2) 
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• Procedure C (C for center): we increase the size of a given sample 
by adding new sites at the center of the sample. So as L increases, 
the Pi close to the boundaries remain unchanged and sites with new 
Pi's are introduced at the center, i.e. further and further from the 
boundaries 

• Procedure B (B for boundary): we add the new sites at the two 
boundaries, whereas the center of the sample remains unchanged. 

For each sample of size L, we make Monte Carlo simulations for different 
values of a at fixed (3 — and different values of a. Figure ^ shows the 
steady state density (■^) where N is the total number of particles in the 
system as a function of a. The density is averaged over typicaly 3 x 10^ 
updates per site. Figure 
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corresponds to three different samples when 



the size is increased according to procedure (C); figure 1(b) to a single 
sample when the size is increased by the sides (procedure (B)), and figure 
1(c) shows the pure case {pi — \). As a varies, one sees a transition from 
a low density for small a to a high density for large a. The slope at the 
crossing point becomes steeper and steeper when the size increases; the 
position of the value of a where the slope is maximum is sample dependent 



but size independent in figure 1(a) when the size is increased by the center. 
On the contrary it is size dependent in figure 1 (b) when the size is increased 
by the sides. 

In the same way, we show in figures (O the variance {{N"^) — {N)'^)/L 
as the size is increased. In the pure case (figure 2(c) I, the location of the 
maximum of the curve converges to Oic — f3, the phase transition point of 
the infinite system as the size increases. When the size of the disordered 



samples is increased by the center (figure 2(a) I, the location of the max- 
imum of the curve converges in a way similar to the pure case but to a 
critical value for a which depends on the initial sample. This suggests that 
in the thermodynamic limit, the system undergoes a phase transition for 
a critical value etc which is sample dependent and depends on the hopping 
rates near boundaries. Moreover, as the size dependence of the maximum 



of the curve is the same in figure 1(a) and 1(c) one expects that the tran- 
sition remains first order in presence of disorder. 



When the size of the sample is increased by the sides (figure 2(b) I, the 
location of the maximum does not converge in the thermodynamic limit. 

Figure|2|shows the distribution of the pseudo transition point (the value 
of a for which the variance in figure [21 is maximum) for quenched systems 
of size L = 50 and L = 100. We generated 40000 samples for each size and 
each sample was simulated over 300000 steps per site, after a sufficiently 
long transient time to reach the steady state. The histogram is then ob- 
tained by discretizing the ac-axe in boxes of length 4 x 10^^. In order 
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(c) 

Figure 1: Mean density versus a for a sample of size 11, 21, 41, 81 
and 161. In figure [l^a)] the system size is increased by adding new sites at 
the center (procedure C), and 3 samples are shown. In figure [r(b)j the size 
is increased by adding new sites on the sides (procedure B) and a single 
sample is shown. In figure |l(c)| the pure case (all the pi = 1) is shown 
for comparison. The dotted line shows on the three figures the transition 
point of the pure TASEP. 
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Figure 2: Variance {{N'^) — {N)'^)/L versus a for a single sample of size 
11, 21, 41, 81 and 161. In figure |2(a)| the system size is increased by 
the center (procedure C). In figure 2(b)| the size is increased by the sides 
(procedure B). In figure |2(c)] the pure case (all the pi = 1) is shown for 
comparison. The dotted line shows the position of the transition point for 
the pure TASEP. 



6 



P(a. 

160 
140 
120 
100 
80 
60 
40 
20 



pure system of size 100 (scaled by 0.-5) 
QASEP of size 50 - 



QASEP of size 100 - 




Figure 3: Distribution P{oic) of the pseudo transition point ac for /3 = 0.15 
when the hopping rates pi are given by Q (full line: size L = 50; dashed 
line: size L — 100). The distribution does not shrink as the size increases. 
The distribution (divided by 2) of the pseudo transition point for the pure 
system measured via the same method is shown for comparison. 



to estimate the numerical noise on the determination of the critical-point 
coming from the finite time length (3 x 10^) of our simulation, we include 
the distribution obtained over 20000 samples of the pure system (the dis- 
tribution for the pure case has been divided by 2 to make the top of the 
peak visible in figure |3Jl. 

The distribution in presence of disorder remains unchanged as the sys- 
tem size L increases by a factor 2, which suggests that the same distribution 
would be observed in the thermodynamic limit. As the distribution P{ac) 
in presence of disorder is much broader than in the pure case, its width 
is not due to the limited time of the measurement. Still, all the details 
of the distribution P{ac) cannot be resolved on a scale smaller than the 
width of the P(ac) of the pure system. Therefore longer simulations should 
reveal more structure at smaller scales of the distribution P{ac) (while the 
distribution for the pure case should become narrower). 

Figure ^ tries to evaluate the correlation between the value of ac and 
the hopping rates close to the boundaries. We note by [.J^^ the average 
over the samples that have a critical point ac (up to the discretization of 
the ac-axis). On figures we show the average of some hopping rates 
LPiJoci namely the first pi, the second p2 and the last pl-i, as a function 
of the value of the critical point ac- One sees that the three maxima of the 
distribution of the critical points can be characterized by the value of the 
first hopping rate pi and pL-i- the left peak corresponds to pi = Pmax, 
PL-i = Pmin, the right peak to pi = Pmin, Pl-i = Pmax- The curve of the 
average of the second hopping rate [p2Jqc as a function of ac supports the 
idea that the value oi p2 has a weaker influence on ac than pi or pl-i- 
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Figure 4: [pija^, [p2jac ^^'^ [Pi-iJoc versus the position of the critical 
point Oc- The full curve is the same distribution of the critical point ac as 
in figure 01 We see that all samples which have a pseudo transition point 
Uc in the left peak have pi = Pmax and pL-i = Pmm and all the sample 
which have a pseudo transition point Uc in the right peak have pi — Pmin 

and PL -1 = Pmax- 



Therefore the picture which emerges is that ac is sample dependent and 
the influence of the p^'s on ac decreases with their distance to the boundary 
(which could lead to a fractal P{ac) (smoothed out in our figure El due to 
the finite time length of our simulation)). The fact that ac is strongly 
correlated to the pi close to the boundaries explains why ac converges 
when procedure C is used as in figure 1(a) and why it wanders around 
when procedure B is used. 



From a technical point of view, let us mention that, for each sample, 
we made the measurement for all the values of a in a single simulation. To 
do so, we introduced K different classes of particles [321 ESI j indexed from 
1 to K, and whose dynamics is the same as in the TASEP, except that the 
particles of a higher class K behave like holes in front of particles of lower 
class K' (when K' < K). More precisely, we still impose the exclusion 
rule, so there is no more than one particle on each site of the system; and 
at each infinitesimal time dt, a particle on site i attempts to jump to its 
right neighbor with probability pidt. The jump succeeds only if the target 
site « + 1 is empty or contains a particle of a class higher than the particle 
which attempts the jump. If the jump succeeds, the particles on site i and 
i + 1 are exchanged. 

To simulate K values ai < a2 < ■ • • < aK of a the rule is to introduce 
particles of class k at rate ak — a^^i (with a^ — 0) at the left boundary 
provided that the first site is empty or occupied by a particle of higher 
class. At the right boundary, a particle present at site L is removed at rate 
(3 irrespective of its class. 
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To do our measurement of the number of particles for a given choice of ak, 
we count as particles all particles of class < k and as holes all particles of 
class > fc + 1 . 

In this work, we have given numerical evidence that for the ASEP with 
open boundaries, the location of the first order phase transition between 
the low and high density phases becomes sample-dependent in presence 
of quenched sitewise disorder. This seems to be a property specific to 
non-equilibrium systems and probably to the one dimensional case as the 
free energy of equilibrium systems (with short range interactions) is self 
averaging and the distribution of the pseudo critical point shrinks as the 
system size increases |34) . 

As in the pure case the transition is controled by the boundary param- 
eters a and f3, our observation that the location of the phase transition is 
influenced by the hopping rates on the sites close to the boundaries is not 
a surprise. In higher dimension, there is a large number of sites closest to 
the boundaries and so their effect should be averaged out. 

It would be interesting to see how the other properties of the phase 
diagram of the TASEP would be affected by site disorder: how the second 
order transitions are modified by the disorder, how the domain wall (uni- 
formly distributed over the sample in the pure case is distributed in 
the disordered case when one sits along at the first order transition point. 

It would be also interesting to see whether these sample dependent 
phase transitions could be seen in other non-equilibrium one dimensional 
systems OZ] when disorder is introduced. 
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